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ABSTRACT 

For over a decade, it has been debated whether the concordance ACDM model is consistent 
with the observed abundance of giant arcs in clusters. While previous theoretical studies have 
focused on properties of the lens and source populations, as well as cosmological effects such 
as dark energy, the impact of initial conditions on the giant-arc abundance is relatively unex- 
plored. Here, we quantify the impact of non-Gaussian initial conditions with the local bispec- 
trum shape on the predicted frequency of giant arcs. Using a path-integral formulation of the 
excursion set formalism, we extend a semi-analytic model for calculating halo concentrations 
to the case of primordial non-Gaussianity, which may be useful for applications outside of 
this work. We find that massive halos tend to collapse earlier in models with positive /nl, 
relative to the Gaussian case, leading to enhanced concentration parameters. The converse 
is true for /nl < 0. In addition to these effects, which change the lensing cross sections, 
non-Gaussianity also modifies the abundance of supercritical clusters available for lensing. 
These combined effects work together to either enhance (/nl > 0) or suppress (/nl < 0) the 
probability of giant-arc formation. Using the best value and 95% confidence levels currently 
available from the Wilkinson Microwave Anisotropy Probe, we find that the giant-arc optical 
depth for sources at 2 is enhanced by ~ 20% and ^ 45% for /nl = 32 and 74 respec- 
tively. In contrast, we calculate a suppression of ~ 5% for /nl = —10. These differences 
translate to similar relative changes in the predicted all-sky number of giant arcs. 
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1 INTRODUCTION 

The formation of giant arcs by strong gravitational lensing is re- 
served for the most massive collapsed structures whose statistical 
properties are sensitive to the expansion history and initial condi- 
tions of the Universe. Since the frequency of giant-arc formation 
depends on the abundance and characteristics of galaxy-clusters 
roughly half-way to the sources, it has long been recognized as a 
potentially rich source of information. 

At the same time, the interplay between cosmological effects, 
cluster physics, and the source population makes their disentan- 
glement non-trivial. The difficulties have been brought to light 
for ov er a decade following the initial claim of iBartelmann et aU 
( Il998l) that ACDM predicted approximately an order of magnitude 
fewer arcs than seen in observations. As subseq uent studies (e.g. 
IZaritskv & Gonzale jl2003l: iGladders et al. 



early observations of lLe Fevre et aLld 19941) . the giant-arc problem. 



l2003h corroborated the 



as it became known, generated considerab le interest because itindi- 
cates one of the following, or both: 1) The lBartelmann et al.l l ll99^ 
analysis was missing a crucial combination of properties exhibited 
by real cluster-lenses and the source population. 2) The concor- 
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dance cosmology is inconsistent with the observed abundance of 
giant arcs. 

Since the first possibility seems most probable, a large amount 
of work has been dedicated towards understanding the most im- 
portant characteristics of arc-producing cluste rs, and how they may 
differ from the general cluster population (e . g . iHennawi et alj2007l : 
iMeneghetti et aLlbOld iFedeU et al.ll2010l) . Other studies focused 
on effects that were not captured in early simulations. For exam- 
ple, it has been shown that artificially populating simulated clusters 
with galaxies in general does not significantly enhance the prob- 
abilit y of giant-arc formation dFlores et al.ll200ol : IMeneghetti et al.l 
I2OOOI) . On the other hand, the mass contribution of central galax- 
ies appears to have a significant effect, th ough not enough to 
entirely resolve the Bartelman n et al] (19981) disagreement alone 
( IMeneghetti etal. 2003; Dalai et alj|2004h . The effects of baryonic 
physics, such as cooling and star formation, on central mass distri- 
butions have also been investigated. The steepening of central mass 
profiles due to baryonic cooling may enhance lensing cross sec- 
tions by factors of a few dPuchwein et al.ll2005l : IWambsganss et al.l 
|2008| : IRozo et al.l2008l. see lMead et alJd201(]|) for a study including 
feedback from Active Galactic Nuclei). 

The properties of background galaxies are equally im- 
portant. The probability of giant-arc formation increases with 



© 0000 RAS 



2 D'Aloisio & Natarajan 



source redshift, making the overall abundance sensitive to uncer- 
tainti es in the high-redshift t a il of the source-redshift distribu- 
tion dWambsganss et alj |2004 lOalal et al . 2004; Li et al. 20oi). 
Moreover, failing to accurately model source sizes and elliptic- 
ities in simulations can alter the expected abundances signifi- 
cantly |Miralda-Escudelll993LlBartelmann et al.|[l995l:IOguril2002l : 



lKeetorJl200ll : IOgurill2002l ; lHo & Whitell2005l ; lGao et aljr2009l) . To 

address these issues, real galaxy images from the Hubble Ultra 
Deep Field have been lensed in recent simulations jHoresh et alJ 
l200il201lh . 

Despite such extensiv e efforts, the status of the giant-ar c prob- 
lem still reniains u n clear (iDalal et all l2004 lLret all I2005L |2006|; 
iHoresh et alj l2005l 120111: iMeneghetti et alj 1201 ih . As iLi et alJ 
( l2006l) and lFedeli et alj | |2008|) point out, the normalization of the 
linear matter power spectrum will play a critical role in determining 
whether there is a giant-arc problem or not. Observations seem to 
be converging on as , « 0.8 jpu et aljlioosi ; Ivikhlinin et al.ll2009l ; 
iKomatsu et alj201lh . while most numerical studies on the giant-arc 
abundance to date have assumed cts = 0.9. It is likely that adjust- 
ing (Tg from 0.9 to 0.8 will lower the predicted giant -arc abundance 
significantly, incr easing tension with observations dLi et al.ll2006l : 
iFedeli et al.ll2008l) . 

With the above caveat in mind, it is still possible that the 
cosmological model may have at least a partial role to play. In 
argu ing that the giant-a rc problem may be unavoidable if ag ~ 
0.8, Ipedeli et"aLl | |2008|) mention in passing that early dark en- 
ergy or non-Gaussian initial conditions may provide "a way out." 
Should such a scenario present itself, the effects of dark en- 
ergy on giant-arc statistics have been well investigated in the past 



le past 



(Bartelmann et al. 2003:' Macci6l2005l : lMeneghetti et alj2005a 

[pedeli & Bartelmann 20071) . On the other hand, the possible effects 
of non-Gaussian initial conditions have not been quantified to date, 
which is the main motivation for this paper. 

We expect primordial non-Gaussianity (PNG) to affect the 
probability of giant-arc formation in at least two ways. Pirst, PNG 
can lead to an enhanced or diminished abu ndance of galaxy clus- 
ters, depending on th e particular model ( e.g. lMatarrese et alj2000l ; 
ILo Verde et ail I2OO8I : lOalal et al.l r2008h . which would lead to a 
change in the number of supercritical lenses that are available 
in the appropriate redshift range. Secon dly, PNG is expected to 
influence the central densities of halo s ( Ayila-Reese et al.ll2003l : 
lOguri & Blandfordll2009l: [Smith et aLlUoiOl) . Since lensing cross 
sections are sensitive to central densities, we expect correspond- 
ing changes in them as well. If a cluster-lens cannot produce arcs 
with length-to-width ratios above some threshold, then its cross 
section for giant-arc production is zero. Roughly speaking, this cor- 
responds to a minimum mass required to produce giant arcs. Owing 
to the effects on central densities, we expect PNG to alter this min- 
imum mass threshold as well. 

A secondary motivation for this work is the question of 
whether giant-arc statistics can potentially serve as a small-scale 
observational probe of PNG. The statistics of rare collapsed struc- 
tures are particularly sensitive to the nature of the primordial den- 
sity fluctuations. Giant arcs are even rarer events and their occur- 
rence is sensitive to subtle changes in the properties of lenses. One 
might therefore expect the effects of PNG to be somewhat ampli- 
fied. Oguri & Blandford (2009) suggest that the statistics of lenses 
with large Einstein radii may be a useful probe of PNG. Here, we 
continue their line of investigation by considering how PNG influ- 
ences giant-arc abundances. We note that the prospect of using arc 
statistics to constrain PNG must be tempered by the considerable 
uncertainties described above. 



The remainder of this work is organized as follows. In ij2]we 
briefly summarize the excursion set formalism and its path integral 
extension for non-Gaussian initial conditions. In ijs] we present a 
semi-analytic calculation quantifying the impact of PNG on the in- 
ner densities of halos. We also compare our calculation to some 
recent simulation results. In Sj?] we calculate the corresponding im- 
pact on the cross section and minimum mass for giant-arc produc- 
tion. In ijs] we present the main results of this paper. We calculate 
changes in the giant-arc optical depth due to PNG. Finally, we offer 
concluding remarks in ij6] 

In what follows, we assume a ACDM cosmology with pa- 
rameters fi™ = 0.27, = 0.73, fib = 0.046, Hq = 
lOO/i kms~^ Mpc"^ (with h = 0.7), Us = 0.97 and ag = 
0.82, consistent with seven- year Wilkinson Micro wave Anisotropy 
Probe (WMAP) constraints ("Komatsu et alJ jgOl ll) . We use the lin- 
ear power spectrum of (.Eisenstein & Hu.l999 



2 THE EXCURSION SET FORMALISM WITH PNG 
2.1 The Gaussian and Markovian case 

In this section, we briefly summarize the excursion set formal- 
ism originally developed for the case of Gaussian initial conditions 
(Press & Schechter 1974; Bond et al. 1991) . For more deta ils, we 
refer the reader to the pedagogical review of lZentneri(l2007l) . 

At its root, the excursion set formalism is a model for es- 
timating the statistical properties of the density field (including 
non-linear growth) using the linearly extrapolated field. The cen- 
tral quantities in the formalism are the density contrast, smoothed 
about some point x, 



<5(x,i?) 



/.V 



,R)5{y 



(1) 



where Wf is a filter function with smoothing scale R, and the un- 
smoothed density contrast is 5(x) = (p(x) — p)/p, where p(x) is 
the mass density and p is the mean cosmic density. In what follows, 
we will exclusively deal with the smoothed density contrast. From 
here on we suppress the x and R dependence for brevity, with the 
understanding that we mean the smoothed quantity. 

The density field at some early epoch is linearly extrapolated 
to a later epoch^ . Working in Lagrangian coordinates, the density 
contrast is smoothed at some large scale around a fiducial particle 
and the variance is calculated. 



S{R) 



\R) = 



(2^ 



P{k)W]{k,R), 



(2) 



where P{k) is the linear matter power spectrum. A useful choice 
for Wf is the coordinate-space top-hat filter, with Fourier Trans- 
form 



Wf{k,R) = 3 



sin(fci?) — kRcos{kR) 

{kWf ■ 



(3) 



The scale of the filter function is decreased and the cor- 
responding 5 and S are again calculated. This procedure is re- 
peated many times, foiming a "trajectory" in (5, 5)-space. When 



^ It is convenient and almost universal in the literature to extrapolate to 
the present day. In what follows, we shall adopt this convention, where the 
redshift dependence is absorbed into the linear threshold for collapse, so 
that (5c (z) &c{z)/D{z) 
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the smoothed density contrast first exceeds some threshold 5c, set 
by a physical model for collapse, the fiducial particle is assumed to 
reside within a halo with mass set by the filter scale R. 

The rate dF/dS that trajectories first cross the barrier in the 
interval S and S + dS is assumed to be equal to the fraction of 
mass contained within halos in the coffesponding mass interval. 
The mass function may therefore be obtained from 



that Sn — en. The probability density in the space of trajectories 
may be written as 

WiSo; Si,..., 5„; S„) = (SoiSiSi) -Si)... fc(5(S,0 - S,,)) 

(8) 

where Sd is the Dirac delta function. The integral representation of 
the Dirac delta function, 



dn{m, ; 
dM 



-dM 



£_ 

AI 



dF 
dM 



dM. 



(4) 



In the excursion set model, the problem of calculating the halo 
mass spectrum is equivalent to determining the distribution of first- 
crossing scales. 

In the specific case of Gaussian initial fluctuations and the 
sharp k-space filter, where Wf{k, R) — 6(1 — k/R), the task is 
simplified considerably and is equivalent to the clas sic problem of a 
Mark ovian random walk with an absorbing barrier dchandrasekhaj 
Il943h . In this case, Il{S, S)dS, the probability density for a tra- 
jectory to obtain a value between S and 5 + d(5 at S, satisfies the 
Fokker-Planck equation. 



'dS 



2 952 ' 



(5) 



with boundary condition Il{Sc, S) — 0. The cumulative probability 
is then given by 



FiS,S) 



n((5,5) dS, 



(6) 



from which the first crossing rate may be obtained by differentia- 
tion. However, one problem with the use of (|5) is that mass asso- 
ciated with the sharp k-space filter is not well defined. A common 
procedure, but one that is no longer necessary (see the next section 
and references therein), is to use the sharp k-space filter in deriva- 
tions and at the end replace it with the form of the coordinate-space 
top-hat. 

A quantity of particular interest for this work is the conditional 
probability that a trajectory will first cross the barrier 5c in the finite 
interval Si to ^2 after having passed through the point (Si , (5i ) . In 
the Gaussian and Markovian case, this cumulative probability is 
given by 



F{S2\Si,Si) =erfc 



Sc — Si 



v/2(S2-5i) 



(7) 



Since equation (O has been used to define the collapse redshift in 
empirical models for halo concentration values, it will serve as the 
starting point for our investigation into how concentration values 
are modified in the case of PNG. 



Sd{5) 



dX -ixs 
2^' ' 



is then used to write 



W{5o;5i,...,Sr,;Sn) = f 



(9) 



(10) 



■y" ,.-.T" , Aj, ...Aj (Si, ...Si )a 



where the bracketed quantities are the connected correlators of the 
smoothed density field, and we have used the notation: 



VX : 



dXi 

2-K 



dX„ 
2tv ■ 



(11) 



The discretized version of the probability density 11 is given by 

rSc 



Tlc{5o; Sn, Sn) 



dSi . .. dSn-iW{So; Si,. ..,5^; Sn). 



(12) 

In practice, F or dF/dS is computed directly by plugging 
into ^ and taking the limit as e — >■ 0. One advantage of for- 
mulating the model in the above way is that non-Markovian effects 
that arise in the case of t he coordinate-space top-ha t filter can be 
treated perturbatively (see lMaggiore & Riottoll20I0ah . A more im- 
portant advantage from our perspective is that the formalism is no 
longer limited to the case in which the higher-order connected cor- 
relators vanish. Models of PNG, which are characterized by higher- 
order connected correlators, may therefore be treated in a fully self- 
consistent way. 



2.3 The local model of PNG 

A common way to parametrize PNG is through the addition of a 
quadratic term in Bardeen's gauge-invariant potential to the usual 
Gaussian piece. 



<E> = "I>G + /nl * [^G - {^g}] 



(13) 



where $g corresponds to Gaussian perturbations and * denotes a 
convolution. The main consequence in assuming perturbations of 
the form ([13} is a non-zero 3-pt correlation function, 



2.2 Generalization to non-Gaussian initial conditions 

We now briefly summarize a pat h integral for mulation of the ex- 
cursion set model, developed bv iMaggiore & Riotto (2010a, b e), 
which has the advantage that it can be applied to a more general 
set of initial conditions, as well as filter functions. We utilize their 
formulation in the next section to estimate the effects of PNG on 
halo density profiles. 

The starting point is to discretize the "time" interval [0, S] so 



($(ki)$(k2)$(k3)) = (27r)'fc(ki + k2 + k3)B*, (14) 

where 5$ = _B#(fci, fe, fcs) is the primordial bispectrum. For sim- 
plicity we consider only the "local" model, where /nl is a constant 
parameter and, to a good approximation, the bispectrum takes the 
shape 



B^{ki,k2,ki) = 2/NL[P#(fcl)P*(fc2) H-P*(fcl)P*(fc3) 

+P#(fc2)P<|.(fc3)]. (15) 
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Here, P# ~ is the power spectrum of the primordial po- 

tential. We reserve investigations of arc statistics with other bis- 
pectrum shapes and scale-dependence for future work. The form of 
the connected three-point coiTelator of the smoothed density field, 
which enters equation dlOt . and its derivative in the local model are 
conveniently summarized in the appendix o f lSmithetal.l l l201(]|) . 

Current constraints on /nl have been obtained from the 
WMAP ye ar seven analysis, whe re the best value is /nl = 32 ± 21 
(68% CL. iKomatsu et alj|201ll) . At the 95% level, /nl is con- 
strained to be —10 < /nl < 74. Comparable constraints on 
PNG have also be en obtained from the Sloan Digital Sky Survey 
jSlosar et al.l2008h . Exploiting the strong impact of PNG on large- 
scale clustering, they obtain —65 < /nl < 70 at the 95% con- 
fidence level. These results are combined with the WMAP con - 
straints to obtain -5 < /nl < 59 (95% CL. lKomatsu et alj20l"ll) . 



3 HALO DENSITY PROFILES 
3.1 The Navarro-Frenk-White model 

The average radial density profiles of halos are well described by 
the profile: 



p{r) 



(16) 



where ps and are strongl y correlated parameters co rresponding 
to density and radius scales jNavarro et al .11 19961 ll997L NFW from 
here on). A convenient way to characterize the density profile of a 
halo is through its concentration parameter, C200 defined as C200 = 
r20o/ra- Unless otherwise stated, we use the convention that the 
mass of a halo is defined by A/ = 200pc(2)47rr2oo/3, where pc(-z) 
is the critical density of the Universe. 

Several algorithms have been developed for obtaining typical 
concentration values for a halo of a given mass and redshift. These 
algorithms generally make use of the apparent connection between 
the central density of a halo and the mean density at t he time of for- 
mation . The first such prescription was developed bv lNavarro et al.l 
( Il997l) . They define the halo collapse redshift Zc as the time at 
which a fraction Fc of the final halo mass M is contained within 
progenitors that are at least as massive as some smaller fraction / 
of M. In their definition, the value of Zc is given implicitly by the 
expression. 



erfc 



5c{Zc) — Sc{z) 

sj2 [S{fM) - S{M)] 



(17) 



which is readily obtained employing the standard Press-Schechter 
formalism. The characteristic density of the halo is assumed to be 
proportional to the mean density at Zc, 



Ss = ps/pc{z) = C f^{z) 



1 + Z 



(18) 



where C is a constant of proportionality that is calibrated by simu- 
lations. NFW suggest using Fc = 1/2, / = 0.01 and C = 3 x 10^. 

Although successful at 2: = 0, subsequent studies showed that 
the NFW prescription above over-predicts con centrations at higher 
redshifts jBullock et aDl200ll : lEke et al.ll200lh for all mass scales, 
and particular ly so for galaxy cluste rs. A lternative algorithms were 
developed by iBullock et ZlfeOOlh and lEke et al.l ( [200 li) to more 



accurately capture t he steeper scalin g of concentration with red- 
shift. More recently. lOao et al] ( l2008h found that adjusting the val- 
ues of Fc =0.1 and C = 600 significantly improves agreement 
with the redshift evolution of cluster conc entrations found in sim- 
ulations^. In what follows, we utilize this lCao et al.l ( I2OO8I) modi- 
fication to the original NFW prescription. The main motivation for 
employing this approach is that it is straightforward to calculate the 
analogous implicit equation for Zc in the case of PNG. 



3.2 The impact of PNG on halo density profiles 

Cosmological simulations in dicate that PNG has an impact on 
the d ensity profiles o f halos jAvila-Reese et al. I I2OO3I J Smith et al.l 
I2OIOI) . Most recently, ISmith et alj ( I2OI0I) found substantial differ- 
ences relative to the Gaussian case in the central regions of ensem- 
ble averaged density profiles. Since strong lensing is very sensitive 
to the central densities of halos, these effects can have a significant 
impact on giant-arc production. To supplement the limited numeri- 
cal analyses available to date, we perform a semi-analytic calcula- 
tion to better quantify this effect. 

Our approach relies on the connection between halo concen- 
trations and their formation redshifts. Using the techniques summa- 
rized in Section|2j2l we quantify how halo formation times change 
for the case of non-Gaussian initial conditions. Interpreted from 
the Gaussian and Markovian excursion set point of view, equa- 
tion ( I17t represents the cumulative probability distribution for the 
ti me at which a single tr ajectory was a fraction / of its final mass. 
As lLacev& Cold d 19931) point out, this does not, strictly speaking, 
yield the distribution of halo formation times. Rather, it represents 
a single progenitor which may or may not correspond to the main 
parent halo in the merger history. However, in what follows, we are 
not concerned with this technical detail. First, it is non-trivial to 
satisfactorily define halo formation times within a simple analytic 
formalism. Second, our primary aim is to estimate differences due 
to PNG relative to the Gaussian case. Since the NFW prescription 
(with slightly modified parameters), based on equation ( 117) , has 
been reasonably successful in its agreement with simulation results, 
particularly in the mass range of interest for this work, it suffices to 
determine how equation ( |17t changes under the influence of PNG. 

A key difference between the Gaussian (and Markovian) and 
PNG cases is that, in the latter, the probabilities for a trajectory to 
propagate from the origin to a point {Sm , <5m) an d from (Sm, 5m ) 
to (Sn,S„) are no longer independent. FoUowing lMa et al. 
the discrete expression for the conditional probability of interest is 
given by: 



/^^ dSi . dSm-i /^^^ d(5m+i . . . dSn-iW {Sq; Si,..., Sn\ Sn) 



W{So; Si, ■ ■ ■ , Sm', Sm) 



(19) 



The above exp ression differs from the conditional probability in 
iMa etalJ ( l2010l) since we incorporate two different barriers through 
the integration limits, Sm = Sc{z2) and Sci = Sc{z\), where z\ > 

Z2- 



2 iDuffv et alj i2008l) sugge s t usin g C = 200. However we find that C = 
600, as found bv lGao etalj <2O08l) . better matches the lOuffv etZI j20O8h 
simulation results - even with the lower value of eg = 0.796. It appears 
that the scaling of concentration with erg is well captured by the modified 
NFW prescription with = 0.1 and C = 600. 
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Figure 1. The effect of primordial non-Gaussianity on halo concentration. 
Top panel: as a function of mass for a fixed redshift of ^ = 0.4. The top 
and bottom dashed lines con'espond to /nl = +100 and —100 respec- 
tively. For comparison, the dotted lines show a si mpler calculation based on 
lOguri & BlandfordI j2009l) (but modified by the iGao et all bOOSh parame- 
ters) which generally agrees well with our results. Bottom panel: same as 
the top panel but as a function of redshift for a fixed mass of 10^^ Mq . 



The cumulative probability of interest is 

F{Sn\5m,Sn,) ^ 1 - f ^ d5„ P {5,, , S„\5„^ , Srr^) ■ (20) 

Note that the Gaussian and Markovian result can be recovered us- 
ing the factorization property, 

^^^■"(^o; 5i, . . . , 5„; Sn) = W'^iSo; <5i, . . . , S„) 
xW^'^iS„,-Sm+i,...,5„;Sn). (21) 

Here, the "gm" superscripts stand for Gaussian and Markovian. 
Plugging equation l l21b into l ll9t . and setting equation l l20t equa- 
tion to Fc yields (TT}. 

The evaluation of equations l |19t and l i20t for the case of PNG 



is quite technical. We provide details in the Appendix. For sim- 
plicity, we do not consider non-Markovian corrections due to the 
coordinate-space top-hat filter in this work. The end result is given 
by equation ( IA35b . which has the form of the Gaussian and Marko- 
vian expression plus terms involving the connected correlators of 
the linear density field, smoothed at the two different mass scales. 
For the case with PNG, equation ( IA35t is set equal to i^c = 0.1 
and solved implicitly for the collapse redshift. The concentration is 
then calculated from the NFW procedure. 

Figure [T] compares halo concentration values of the PNG 
case with that of the Gaussian case. The dashed curves show 
the results from our approach. For comparison, we dis p lay a n 
estimate based on the approach in lOguri & BlandfordI (|2009|) . 
They use an effective rescaling of the collapse barrier, 5c{z) — >■ 
Sc{z)y^l — S36c{z)/3, where 53 is the skewness, which is mo - 
tivated by the ma ss function deriv ation of iMatarrese et al. IdlooO). 
We adjust for the lGao et al.l ( 1200 Sl) modifications and neglect halo 
triaxiality. Their approach generally agrees well with our results. 
However, in our calculation, larger mass halos are more affected 
by PNG. Similarly, higher redshift halos display a larger deviation 
from the Gaussian case. 

In Fig u re [2] we compare our calculation to the results of 
ISmith et aU ( |2010|) . which are derived from N-body simulations. 
We show the ratio of density profiles as a function of ra dius for two 
differ ent masses in units of M@/h (see Figure 7 of I Smith et al] 
l2010i) . Some c are needs to be exe rcised due to differing definitions 
of halo mass. ISmith et alj ilOld) use Friends-of-Friends masses. 
They point out a reasonable agreement between these masses and 
masses defined with the 200 x p criterion. We therefore convert'^ 
our masses and concentrations to be consistent with the 200 x p 
definition in Figure|2] 

The solid lines in Figure |2] show our calculation while the 
dashed lines correspond to lo g-linear model fits to the ensemble 
averaged density profiles in the lSmithetaLl(l2010l) simulations (see 
Figure 7 of their paper for simulation data). Note that their results 
are obtained from stacking within bins, whereas our results are cal- 
culated for the average masses of the bins. The vertical dotted line 
shows their softening length of 40 kpc. Models with positive /nl 
yield enhanced central densities and vice versa. Our calculation 
agrees reasonably well with the numerical results, particularly in 
the inner regions of halos, which are most important for giant-arc 
production. 



4 THE GIANT-ARC CROSS SECTION 
4.1 Gravitational lensing 

In the case where the deflecting mass distribution is localized rela- 
tive to the cosmological distance scales involved, gravitational lens- 
ing is described to a good approximation by: 



f3 = e-a{e), 



(22) 



where /3 and 6 are angular positions in the source and image planes 
respectively, and a is the deflection angle, which may be obtained 
directly from the lensing potential tp through 



^ A procedure for co nverting between mass definitions can be found in 
ISmith & Watt3 i2005l) 
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r (Mpc/h) 

Figure 2. The effects of primordial non-Gaussianity on halo density pro- 
files. We compare our semi-analytic calculation (solid lines) for two differ- 
ent masses (in this plot, M i s give n in Mq //i) to recent results from N-body 
simulations in lSmith et alj ilOlOl) . This plot may be compared to Figure 7 
of their paper. The dashed lines show their log-linear model fit. The vertical 
dotted line shows their softening length (40 kpc). Models with /nl > 
yield enhanced central densities and vice versa. These changes can i mpac t 
the cross section for giant-arc production. Note that the Ism ith et al. ( 2 0ld) 
results have been obtained from halos stacked within mass bins, whereas 
our results are for the average masses of those bins. 



side of equation ([27} form the tangential critical curve; named as 
such to reflect the typical orientation of nearby images. Conversely, 
the roots associated with the second factor form the radial critical 
curve. 

The source plane locations associated with the critical curves 
form the caustics of the lens mapping. Background galaxies that re- 
side sufficiently close to the caustics may be lensed into highly dis- 
torted images. The area in the source plane corresponding to images 
with length-to-width ratios above some threshold is the cross sec- 
tion for giant-arc production, which we will denote with aa from 
here on. 



4.2 Lensing by pseudo- elliptical NFW halos 

We utilize a pseudo-el liptical NFW lensing potential developed by 
iGolse & Kneihl i2002l) . Below, we summarize the circularly sym- 
metric NFW lens and the procedure of Golse & Kneib (2002) to 
obtain the pseudo-elliptical extension. 

Given that the critical curves are determined by the condition 
l |27| l, the relevant quantities for our purposes are the convergence 
and shear. For a circularly symmetric NFW lens, the former is given 
by 

n{x) = 2Ksf{x), (28) 
where Ks = Pcrs'Sc^, and x = r/rs. The shear is 

7(x) = - f{x)^ , (29) 

where we define 



a = Vip. (23) 
The local distortion of images can be quantified by the Jacobian, 



1 



K — 71 
-72 



1 



-72 

K + 71 



(24) 



Here, the convergence k is related to the lensing potential through 
the two-dimensional Poisson equation. 



= 2k. 



(25) 



The components of the shear, 71 and 72, are also given by second 
derivatives of the lensing potential. 



71 = 2 (^.11 

72 = •'P,12- 



(26) 



In this work we are interested in images that are highly dis- 
torted. Such cases typically occur near the critical curves of the lens 
mapping i22\ , which are formed by the points where the Jacobian 
is singular. These points satisfy det {df3/dG) = 0, or 



1 - 



, ^ arcch- 



^11 



l2-l 



and 



In- 



= arcch- 



1 + lni 



{x<l) 
ix = l) 
ix>l) 



{x<l) 

{x = l) 
{x>l) 



(30) 



(31) 



In the procedure of iGolse & KneibI 1 I2OO2I) , the coordinate transfor- 
mation. 



e Xi 



Xi - 
X2 — > Vl + e X2, 



(32) 



is applied to the above equations in order to introduce ellipticity 
and generalize the lensing potential. The ellipticity parameter e is 
related to the major and minor axes (a and b respectively) of the 
iso-potential ellipses by 



(l-K-|7|)(l-«+|7l) =0, 



(27) 



where I7I is the magnitude of the complex valued shear 7 = 
71 + ^72 ■ The roots associated with the first factor on the left-hand 



62' 



(33) 



As lGolse & Kneihl j2002l ) point out, the contours of the lensing po- 
tential become more "peanut" shaped as the ellipticity is increased 
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to e ~ 0.3 and beyond. This regime of hiigii e values can lead to neg- 
ative mass densities at larger radii. For these reasons, and for agree- 
ment with recent ob servational results (e.g. iLimousin et al.|[20Q8l : 
lNewmanetal]|2009l) . we restrict ourselves to values of e = 0.1 
and 0.2 in the calculations below. 

Under the transformation i32\ , the convergence and shear be- 
come, 



Ke(x) = K,{xt) + ecos(2</)£) 7(a;e) 

[hi^{Xi) - sm{2(l>i)'y'^{xe)] , 
where we use the polar coordinates. 



(34) 



(35) 



x^ = \Jx\^ + xl^ = y^(l - e) + (1 + e) X 
— arctan(a;2£/a::ie). 



(36) 



4.3 The cross section and minimum mass for giant-arc 
production 

Having estimated the impact of PNG on halo density profiles 
in Section 13.21 the primary aim of this section is to explore 
the corresponding changes in giant-arc cross sections. Rather 
than use computationally expensive ray tracing techniques (see 
ISartelmann & Weiss 1 994|, f or example), or the surface integral 
method o f fpedeli et al.l hooeh to calculate cross sections, we use a 
simple approximation used bv lBartelmann et al] ( l2003l) which cap- 
tures the scaling of aa with mass and redshift (which we have in- 
dependently checked). 

Let {±0a, 0) and (0, ±9b) be the locations where the tangen- 
tial critical curve intersects the coordinate axes. We use equations 
( I34t and i35l with the first factor in (I27t t o determine these lo- 
cations. Following iBartelmann et al.l l l2003h . we assume that the 
giant-arc cross section scales approximately with the area enclosed 
by the tangential critical curve, aa ~ OaOb- Since our main goal 
is to calculate changes relative to the Gaussian case (i.e. ratios of 
quantities), we do not need to know the constant of proportionality. 
Hence, we assume that aa ~ OaOb with the understanding that aa 
is not the absolute cross section, but merely an approximation up to 
some multiplicative constant. 

Figure |3] shows the ratio of non-Gaussian to Gaussian cross 
sections as a function of /nl for several different halo masses. For 
/nl 7^ 0, we use the concentration values obtained from equation 
( |A35b , as described in Section[3] We use a lens redshift of zi = 0.4, 
an ellipticity of e = 0.2, and a source redshift of Zs = 1.82, which 
is the median redsh ift observed in the Sloan Giant Arcs Survey 
teavliss et al .1201 Ih . As we showed in Section [X2l halos in models 
with /nl > have enhanced central densities. Such clusters there- 
fore have greater cross sections for giant-arc production. The con- 
verse is true for /nl < 0. Figure[3]shows that PNG can enhance (or 
decrease) giant-arc cross sections by up to 20% for |/nl| ~ 100. 
For the large masses considered in this work, the effects are only 
mildly dependent on the mass due to the fact that, in those cases, 
giant arcs tend be located at larger distances from the center. The 
relative changes in the central densities therefore have less of an 
impact as M is increased. 

In practice, arcs are only considered "giant" if their length-to- 
width ratio exceeds some threshold value. Therefore, if a lens is not 
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Figure 3. The ratio of giant-arc cross sections in the case of non-Gaussian 
and Gaussian initial conditions. Halos have enhanced central densities in 
models with /nl > 0. Their giant-aix cross sections are therefore increased 
relative to the Gaussian case and vice versa. 



capable of producing arcs a bove the given threshold , then its cross 
section is taken to be zero. IBartelmann et al.l | |2003|) incorporate a 
condition which is meant to emulate this threshold. If the major 
axis of the tangential critical curve is below ^min, then the cross 
section for giant arcs is set equal to zero. Following their model, 
we use a fiducial value of ^min = 10", but explore the impact of 
using Smin = 5" ou the main results of this paper in Section ISTI 

For given lens and source redshifts, 6*111111 translates to a mini- 
mum mass, Afniin, for giant-arc production. In Figure |4] we show 
Mniin corresponding to Omin = 10" as a function of lens redshift. 
Since halos have a higher central density for /nl > relative 
to the Gaussian case, halos can be less massive and still meet the 
^min = 10" condition. The converse is true for /nl < 0. 



5 RESULTS 

5.1 Giant-arc probabilites 

The probability for a source at redshift Zs to produce giant arcs is 
given by the optical depth**, 



az 



AM — a^{M,z), (37) 



where dV/dz is the comoving volume element. An /AM is the halo 
mass function, and Mmi-a is the minimum mass to produce giant 
arcs (see Section l43t . 

In the previous sections we have discussed the effects of PNG 
on the cross section and minimum mass. The final ingredient for 
our simple model is An /AM. W e utilize one of the mass func- 
tions tested against simulations in ISmith et alj bOld) . which is a 



Note that we have utilized the analytic approximation, a a ~ daOi,- In 
this case, the cross section is in angular units. Note that the angular diameter 
distance to Zs does not appear in equation )37t . 
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Figure 4. The minim um mass for giant-arc production. Following 
iBartelmann et al.l )2003l) . we impose a threshold for what is considered a 
"giant" arc that is based on the major axis of the tangential critical curve. A 
given lens has a non-zero cross section for giant-arc production if the major 
axis is above f^min ■ 

The ^?min threshold effectively results in a minimum 
mass, Afjnin. Models with positive /nl lower Afniin due to the enhanced 
central densities and vice versa. The changes in M^i^ play an important 
role when integrating over the lens population to obtain the giant-arc opti- 
cal depth. 



slight modification of the lLo Verde et all (l2008h form ^. Note that 
we must convert masses due to the fact that the M which appears 
in the mass function is defined in terms of 200p rather than 200pc 
(which is the convention that we have used). 

Figure |5] shows the integrand of equation ( I37t . which in our 
calculation is the differential optical depth At /Az up to a multi- 
plicative constant (see Section |43t . For lenses at zi ~ 0.4, PNG 
with /nl = 100 increases the differential optical depth by ~ 50% 
relative to the Gaussian case and vice versa. 

The giant-arc optical depth for the median redshift in the Sloan 
Giant Arcs Survey is shown in Figure |6] We show results in the 
range —100 ^ /nl ^ 100 and for various combinations of e and 
^min- We note that the deviations in t from the Gaussian case are 
due to the combined effects of modified central densities and halo 
abundance. For example, in the case with /nl > 0, central den- 
sities are enhanced and the abundance of high-mass halos is in- 
creased, which can boost the giant-arc optical depth substantially. 
For the bes t estimate of /nl = 32, obtained from the WMAP year 
7 analysis jKomatsu et alj|201 ih , we estimate a ^ 20% increase in 
the giant-arc optical depth relative to the case with /nl = 0, for 
sources at Zs = 1.82. Using the WMAP 95% confidence levels, we 
calculate a —5% and +45% change for /nl = — 10 and /nl = 74 
respectively. PNG makes less of an impact for lower values of 0min ■ 
This is due to the fact that lower ^min corresponds to lower Afmin- 
In this case, low-mass halos whose abundance and central densities 



^ The expression obtained with the Edgeworth expansion in lLo Verde et alj 
)2008l) , which is used as a correction factor to the mass function, may also 
be obtained with the formalism described in Section 12. 2| It is derived from 
leading order terms with the sharp k-space filter (see llVIaggiore & Riottol 
l2010ch . 



Figure 5. The differential giant-arc optical depth (up to a multiplicative 
constant) as a function of lens redshift. For lenses at typical redshifts (z; 
0.4), PNG with /nl = 100 enhances the differential optical depth by ~ 50 
per cent relative to the Gaussian case and vice versa 
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Figure 6. The ratio of non-Gaussian to Gaussian giant-arc optical depths. 
We show results for a few different combinations of the lens potential el- 
lipticity, e, and minimum value for the major axis of the tangential critical 
curve, 6*min. 



are less affected by PNG contribute more to the optical depth. On 
the other hand, a higher value of e leads to less relative change from 
PNG. This is due to the fact that the giant-arc cross section grows 
with e, making the relative contribution from central densities less 
important. 

Figure |7] shows the ratio of non-Gaussian to Gaussian optical 
depths as a function of Zs . The effect of PNG decreases mildly with 
the source redshift. At first glance, one might find this surprising; 
PNG has a larger impact on the halo mass function and density 
profiles at higher zi. The overall decrease results from two subtle 
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Figure 7. Same as Figure |6] but as a function of source redshift, and fixed 
values of /nl- Although the optical depth increases with source redshift, 
the impact of PNG (as shown thi'ough the ratio) decreases due to two effects. 
First, the typical distance of giant arcs from the center of lenses increases 
with Zs, making effects on the central densities less important. Secondly, 
for high Zs , the minimum mass for giant-arc production is decreased over a 
wider range. In this case, lower mass halos, whose densities and abundances 
are less influenced by PNG, contribute more to the optical depth. 



Figure 8. The number of giant arcs per unit lens-redshift (up to a mul- 
tiplicative constant). Although our aim is not to calculate the absolute arc 
abundance, the top panel is ai'bitraiily normahzed so that the all-sky num- 
ber of arcs in the Gaussian case (solid line) is 1000, roughly in accordance 
with observational results. The location of the peaks in the top panel illus- 
trate that most arcs will be observed within cluster-lenses at 2; 0.5. The 
bottom panel illustrates that the largest fractional changes due to PNG occur 
in the lower and higher redshifts regimes, where giant arcs are rarest. 



effects: 1) For a fixed mass and zi, the average radius of the criti- 
cal curves (and caustics) grows with Zg. Giant arcs tend to reside at 
larger radii from the cluster core. Hence, the effects of PNG on the 
central densities become less important at larger Zs . 2) The mini- 
mum mass for giant-arc production decreases over a larger range 
of zi when Zs increases. In this case, due to the steepness of the 
mass function, the integral over M in equation J37t receives more 
contribution from dn/dAI in a lower mass regime, where PNG has 
less of an effect. 



5.2 Giant-arc abundances 

Since we cannot calculate the absolute optical depth with the semi- 
analytic approach taken here, predicting giant-arc abundances is 
well beyond the scope of this paper. However, we can estimate the 
relative changes due to PNG which is of great interest. The ex- 
pected number of giant arcs per square degree is 



^ arcs — 



rdzs^T{zs), (38) 
Jo azs 



where dNs/dzs is the differential source density. We use a fixed 
dNs/dzs obtained from the observed galaxy redshift dis tribution 
in the Canada-France-Hawaii Telescope Legacy Survey jFu et al.l 
l2008h . We evaluate the integral ([38} up to Zs = 5. For /nl = 32, 
— 10, and 74, we obtain changes in the predicted number of giant 
arcs per square degree of +17%, —5%, and +41% respectively. 
For an extreme value of value of /nl = —100, the change is 
—43%. In order to get "back-of-the-envelope" estimates of what 
this implies in practice, we use the all-sky extrapolation of ~ 1000 
arcs wit h length-to-width ra t io ^ 10 and R-band magnitudes 
< 21.5 jLe Fevre et ailll994 Isartelmann et ai]|l998l; lOalal et al.l 



l2004 l Fedeli et al. Illooi). Assuming that the theoretical prediction 
for the Gaussian case is of this order, the non-Gaussian cases with 
/nl = 32 and 74 would predict 170 and 410 more giant arcs over 
the entire sky respectively, while /nl ~ —10 would predict 50 
less. 

In figure [S] we plot the number of arcs per unit lens redshift, 
dA^arcs/dz, (again, up to a multiplicative factor). We arbitrarily 
normalize dA'^arcs/dz; so that the all-sky number of arcs for the 
Gaussian case is 1000. From this plot, one may infer changes due 
to PNG in the number of arcs observed in clusters at a given lens 
redshift. The top panel illustrates that most arcs will be observed 
in clusters at zi ~ 0.5. The bottom panel shows that the largest 
fractional changes due to PNG are at the low and high redshift 
ends of the distribution. Both of these effects are mainly due to 
the fact that the minimum mass threshold (see figure |4j increases 
at both lower and higher lens redshifts. In these redshift regimes, 
the lenses containing giant arcs typically correspond to the most 
massive and rarest peaks, whose statistical properties are most af- 
fected by PNG. The largest fractional changes will therefore be in 
the redshift regimes where giant arcs are extremely rare. 



6 SUMMARY AND DISCUSSION 

We have quantified the impact of local primordial non-Gaussianity 
on the statistics of giant arcs in clusters. Our calculations take into 
account changes in both the abundance and central densities of 
clusters due to non-Gaussianity. We quantifi ed the effect on centra l 
densities by extending the a nalytic model of Navarro et al.l d 19971) , 
with parameters modified bv lGao et al.l | |2008|) , for calculating typ- 
ical concentration values as a function of mass and redshift. Our 
approach utilizes a recently developed path integral formation of 
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the excursion set model jMaggiore & Riottoll2010allbll3) . We calcu- 
late corrections to the collapse fraction whic h is used to implicitly 
define halo collapse redshifts in the original iNavarro eFal 
model. We find that massive halos tend to collapse earlier in mod- 
els with positive /nl and vice versa. Since the central densities of 
halos reflect the cosmic density at their collapse epoch, this leads 
to enhanced and suppressed central densities for positive and neg- 
ative /nl respectively. As an example, for a halo of mass lO^^M© 
at z = 0.4, we calculate a ±5% change in concentration parameter 
for /nl ± 100. 

We have compared our estimates of the rel ative change in con- 
centra tion values with a different approachjjy Oguri & Blandford 
( I2OO9I) and recent N-body simulations bv lSmith et al.. (2010.,) . and 
find good general agreement with both. Our approach may be par- 
ticularly helpful in future numerical studies which aim to determine 
how primordial non-Gaussianity impacts halo density profiles. 

The modified halo abundances and changes in central densities 
work in the same direction to either enhance (/nl > 0) or suppress 
(/nl < 0) the probability of giant-arc formation. The effects of 
PNG on halo abundance alters the number of supercritical clusters 
available for lensing. The central densities make an impact in two 
ways: 1) They affect the lensing cross sections. For lenses with 
M ~ lO^^'Mo and zi ~ 0.4, we estimate changes of up to ±20% 
for /nl ± 100. 2) Since more or less mass may be concentrated 
in central regions, the minimum total mass that a lens must have 
to produce giant arcs changes. For /nl ± 100, the effect on the 
mass threshold is or order ±10% for 0.2 < zi < 0.6. Note that 
while our analytic approach requires changes in the mass threshold 
to be imposed by hand, they would appear naturally in ray tracing 
simulations through the altered cross sections. 

We have calculated changes in the giant-arc optical depth rel- 
ative to the Gaussian case. For a source redshift of Zs = 1.82, cor- 
respond ing to the median v alue observed in the Sloan Giant Arcs 
Survey jSavliss et al.l201ll) . the above effects translate to a ~ 20% 
enhancement of the optical depth for the WMAP year seven best 
value of /nl = 32. For the 95% confidence levels of /nl = — 10 
and /nl = 74, we obtain changes of approximately —5% and 
+45% respectively. The relative change in the optical depth due to 
non-Gaussianity decreases mildly with source redshift. This is due 
to the fact that for a fixed lens redshift, the radii of critical curves 
increase with the source redshift, so that the effect on central den- 
sities becomes less important. Also, the minimum mass threshold 
for producing giant arcs is generally lowered as the source redshift 
increases. The optical depth therefore receives more contribution 
from lower mass halos which are less affected by non-Gaussianity. 

We have also calculated changes in the predicted giant-arc 
abundance due to non-Gaussianity. For /nl = 32 and 74, we ob- 
tain 17% and 41% enhancements in the predicted number of gi- 
ant arcs per square degree respectively. In contrast, /nl = — 10 
leads to a 5% decrease. We have examined the number of giant arcs 
per unit lens redshift, finding that most will be observed in cluster 
lenses at zi ~ 0.5, but that the largest fractional changes due to 
non-Gaussianity will be in the low {zi ^ 0.2) and high (zi > 1) 
lens-redshift regimes, where arcs are extremely rare. The statistics 
of giant arcs in clusters solely in those redshift ranges are likely to 
be limited by cosmic variance. 

There are some other ways that non-Gaussianity can influence 
giant-arc statistics, which we have not been able to consider with 
our semi-analytic approach. One possible way is through its influ- 
ence on the clustering of massive halos (e.g. Matarres e & Verd3 
l2008l : iDalal et alj |2008| : ICarbone et af] l2008l : iMcDonall l2008h . 
which can change the role of neighboring structures. Addition- 



ally, since non-Gaussi anity may introduce large-scale corre lations 
between clusters (e.g. iFedeli et al.ll2009l ; ICunha et al.ll2010l) . their 
lensing properties may be more influenced by line-of-sight align- 
ment. These are topics which are best addressed in the future by 
ray tracing through cosmological simulations. 

Another possible way is through halo merger events, which 
have been shown t o be important for a rc statistics in previo us works 
jTorri et al.ll2004l iFedeli et al.ll2006l) . iFedeli et al.l ( |2006|) have de- 
veloped a semi-analytic Monte Carlo method for incorporating 
mergers without computationally expensive ray tracing techniques. 
We did not pursue their approach here due to the fact that non- 
Gaussian initial conditions introduce correlations between scales, 
which complicates the algorithm considerably. However, since non- 
Gaussianity may influence merger rates, one might also expect to 
see corresponding differences in the optical depth relative to the 
Gaussian case. The impact of non-Gaussianity on merger rates 
may b e studied analytically using techniques in iDe Simone et aU 
dioi l|), which we became aware of during the preparation of this 
manuscript. This is a topic that is interesting both in its own re- 
spects and for its implications. 

Non-Gaussianity could also impact the giant-arc abundance 
through additional effects on halo structure and substucture. Our 
zeroth order approximation only considered the impact on halo 
concentrations. For simplicity, we do not consider substructure in 
this work. However, more detailed effects of non-Gaussianity on 
density profiles and the mass spectra of substructure can make a 
significant difference in the lensing properties of clusters. Again 
this a topic that is best treated numerically through high-resolution 
simulations. We note that the above effects would also contribute 
changes to the giant-arc cross section, making changes in central 
densities relatively less important. 

We now discuss the relevance of our results to the giant-arc 
problem summarized in the introduction of this paper. Within the 
restrictive constraints on /nl from WMAP, the effects of local-type 
PNG are relatively modest, as one might suspect, implying that it 
cannot account for an order-of-magnitude giant-arc deficiency, if 
one should exist. On the other hand, recent detailed investigations 
find much less disagreement with observations compared to th e 
original ^artelmann et al. (1998) results (e.g. [Horesh et al.i l201lh . 
We emphasize again that most giant-arc studies to date have used 
(Tg = 0.9, so it is still reasonable to expect some some lev el of dis- 
agreement when adjusted for a lower value of erg = 0.8 \1a et al. 



l2004 IFedeli et alj|2008l) . We examine Figure 10 of IFedeli et a~ 
12008!) to estimate the possible level of disagreement when erg = 
0.9 is adjusted to erg = 0.8. Using the bottom curve corresponding 
to a limiting R-band magnitude of 21.5, and length-to-width ratio 
^ 10, we estimate a ~ 60% decrease in the predicted all-sky num- 
ber of arcs. If this is the case, then local-type non-Gaussianity with 
/nl ~ 32 is on the right order of magnitude (~ 20%) to at least 
help compensate for this deficient. 

We note that additional effects, as discussed above, may con- 
tribute to a larger impact than estimated in this work. Also, other 
models of non-Gaussianity, which correspond to different bispec- 
trum shapes, may result in a larger impact, particularly when 
scale-dependence is introduced. The scale dependence allows non- 
Gaussianity to have a greater influence on smaller scales relevant 
to structure formation, and giant-arc statistics, while minimizing 
effects on s cales relevant for Cos mic Microwave Background mea- 
surements ( ILo Verde et al. ||200i. This is a topic that we are cur- 
rently addressing using the foundations laid in this work. 

Finally, we address the question of whether arc statistics can 
some day serve as a complementary observational probe of PNG, 
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when larger lensing-cluster samples are acquired. Certainly, de- 
tailed numerical simulations would need to be performed in order to 
more precisely quantify the effects of PNG. While the ~ 50% level 
effect for /nl ~ ±100 (local) calculated here seems significant, 
these effects are still relatively weak compared to uncertainties in 
even the more sophisticated numerical methods. Additionally, even 
with detailed simulations in hand, a fai r comparison with obser - 
vations is non-trivial. As pointed out by iMeneghetti et al.l ( 1200 Sh . 
a number of observational effects will need to be accurately sim- 
ulated for the comparison. For example, background noise from 
other photon sources, atmospheric effects and the point-spread- 
function can lead to altered length-to- width ratios and compromise 
arc detectability. In fact, the observed abundance of giant arcs de- 
pends sensitively on the characteristics of the survey. Another dif- 
ficulty stems from the fact that real cluster lenses are selected on 
observables such as X-ray luminosity. It is crucial to accurately 
match the survey selection criteria in simulations since the giant-arc 
abundance can var y significantly depending on the selection limit 
jpedeh et alj20IOl) . Therefore, although giant-arc samples will cer- 
tainly grow in the coming years, there are a number of uncertain- 
ties, both theoretical and observational, that need to be better char- 
acterized and reduced before arc-statistics can serve as a probe of 
PNG. 
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APPENDIX A: THE CONDITIONAL FIRST-CROSSING PROBABILITY 

In this section we outline a derivation of the cumulative, conditional probability i20\ in the case of non-Gaussian initial conditions. We begin 
by considering the conditional probability given by equation \191 . By expanding the exponential in equation ilOl . we can write 



W {60; ■ ■ ■ , S„; Sn) = I 2?Aexp-j iXiSi + ^' XiXjXk {5i5j5k), 



1 " 



where the first and second terms on the right hand side correspond to the Gaussian and first-order non-Gaussian (three point connected 
coTOlator) contributions respectively. The sum in the second term of equation dAlt can be broken up so that 



m — 1 



- ^ {SiSj5k)c didjdk = - - ^ {SiSjSk}c didjdk - - Y^ {SiSjS,n)cdidjd„ 

i,j,k — l i,j,k — l 

1 1 



a3 



(A2) 



m — 1 n 



— l fc— m + 1 
m — 1 n 

z — 1 J — m + 1 



i — 1 j ,k — m-\-l 



(A3) 



1 " 1 " 

- ^ ((5i(5j5fc)c dtdjdk - - ^ (5i5j5m)c didjd„ 

i,j,fc — m + l i.j — m + 1 

1 

i — m+ 1 



(A4) 



To begin, we consider the three terms in 1 IA21 . For brevity, we use the following notation: = W^'^{5m', Sm+i, ■ ■ ■ , 5n; S„). Using the 

factorization property ||211 and the chain rule, we apply the derivatives with respect to 5m to obtain 



^ m— 1 

m— 1 
m — 1 



(A5) 



where we have obtained the first term on the right hand side by combining terms. Plugging this into equation M9\ along with the contributions 
from ( IA3t and ( IA4t . and keeping only terms that are first-order in the connected 3-pt correlators yields 



cll'Jm; 5n\ Sn — Sm) + -Pe ^{5n, Sn\Sm, Sm), 



(A6) 



where we have defined 



^ {^n ) Sn \ Sm •, Sm) 



nf™((5m|(5o; Sm; Sm — So) ' 



(A7) 



with 
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Na 



d6l . . . d(5m-l / d(5m + l . . . d6n-i 
-oo •/ — oo 

m— 1 m— 1 



i,J=l 

m — 1 



m — 1 / r3 \ / £"3 \ 

i=i 



<5m /■<5. 

d5i . . . d5m~i I dS, 

— oc 

m—1 



+I...A5 

n — 1 

^ ^ ,„ ^ vn — 1 n 



fc — m+ 1 



m — 1 n m—1 n 

z — 1 J— m+1 i—l J— m + 1 



<5m r^cl 

d5i...d5m_i / d(5m+i . . . d(5„_i 



i.j — m + l 



i , J — m + 1 

m + 1 i — m+1 

n n 



,2 pif^gm 
m ' ^ m.n 



' m.n 



i — m+ 1 



(A8) 



(A9) 



(AlO) 



At first glance, the evaluation of these expressions appears to be a formidable task. However, we can make progress by usi ng some simplifying 
assump tions involving the connected correlators. Using a procedure similar to the discussion following equation (41) in lMaggiore & Riottd 
( l2010d) . we make the substitutions; 



J2 {s^s,s^) ^ (Si) E 

m— 1 m—1 

^{5Al>^(e)E 

1^1 i^l 
m—1 n m—1 n 

E E ^ (^'rnSr.) E E 

— l fc — m+1 *;j — 1 fc — m+1 

m — 1 n m — ln 

E E {s^w^{s^si)j2 E 

i— 1 j,fc— m + 1 i — l j.fc — m+1 

m. — 1 71 m—ln 

E E {s^s,s„,) ^ {sis^} J2 E 

i—l J— m + l i—l j — m+l 

n n 

E {5.5,6,} ^{SD J2 

i,j,fc— m + 1 i,j,fc — m+1 

n n 

E {5iSjSm) ^ (5m5S) E 



— m + 1 



i, j — m + 1 



E (5i5m} ~ E 
i— m + 1 i — m+] 



(All) 
(A 12) 
(A13) 
(A14) 
(A15) 
(A16) 
(A17) 
(A18) 
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We also make exten sive use of the tricks given by equations (48), (49), and (50) of iMaggiore & Riottol ( 1201 Pel) (see also 
iMaggiore & Riottoll2010al) . However, we will find it necessary to treat cases in which the upper limits of integration also appear explic- 
itly in the integrand. For these cases, it may be shown that 



J2 d5i . . . d6m-i ftM/o^,^ + / d5i . . . d<5^_i dmW§;, 

T J —oo J —oo 



a^nr(5„|5o;<5™;S™-5o) = > / ASx . . . d&m-i d.WI'^ + AS^ . . . dS^-i drr.WI'^ (A19) 



=l2 



i = l 

: J —oc J 1 J—cyo 



= i=l 

+ dSi... dSr^-i 9mVK(f,^, (A20) 

J — OO 

m~-l pSm m — 1 pSjn 

dlni'^{5m\5o;Sm;Sm-So)^ / d5i...d5„,-ididjdkW§^',^ + 3Y^ dSi . . .d5„,^i didjO^W^^^ 

+3 ^ / d5i . . . d5^-i drdl^W^Z + d5i . . . dSm-i d^WiZ- (A21) 

In applying these tricks, we run into the same issue as described in section 3.1 of IMaggiore & Riottol ( l2010ch : namely that some of the 
summations in equations l|A9l > and dAlOl l are up to n and not to n — 1. They are therefore not in the form of equations (48), (49), and (50) in 
IMaggiore & Riottol bOlOdT As they point out, we are ultimately interested in calculating F{S„\Sm., Sm)- This is given by 

F{S„\Sm, Sm) = 1 - d5„ UnScl\5m; 5„- Sn - Sm) - dSn P^^'iSn, S„\5m, Sm) ■ (A22) 

J — oo J — DC 

We will therefore evaluate J^"^ dSn Na{bc) instead of Na(bc) - For terms with the strategy is to substitute the first terms on the 

right-hand sides of equati ons (|A19ll - (IA21b wherever possible. Similarly, for terms involving J]]", we aim to substitute the right-hand sides 
of equations (48) - (50) in lMaggiore & Riottol ( l2010cl) . After a fortunate cancelation of all terms involving integrals over derivatives of VKy ™ 
with respect to Sm, we find that 

dS„ Na = -^dl [nnSm\0,m)]dm [K(m,n)] - dm [nn5m\0,m)]dl [K(m,n)] 



2 



(Si) nr(5„|0,m)a^[(7.(m,n)]. (A23) 



6 

For brevity, we have used the shorthand notation: nf™((5m|0,m) = Ilf^ {Sm\So; Sm', Sm ~ So). We also define Ut{'m,n 
d5„nf™(5ci I^Ti, n). Following a similar procedure for A^;, and Nc yields 

f'' dSr, Nt=- ol [nr(5m|0, m)] del [U,{m, n)] - dm pr (<5m|0, m)] d^i [C/.(m, n)] 

- {SmSn) dm [nf™(5„i|0, m)] dmdci [Ui{m, n)] 

+ ^^|^9ci [U.{m, n)] dSi... dSm-idlWlZ + ''^^dl^ [f/.(m, n)] d5i . . . dSm-idmWl"^ 

^ J -oo ^ J -oo 



+ {S''mS^)dAdm \U.{rn, n)\ I dSl . . . dSm-ldmW§^ (A24) 



I dSi 

■J — oo 



[ dS„ Nc=- -^^2^ nr{Sm\0,m)dldci [U4m,n)] - -^^^ nnSm\0,m)dmd^,^ [U4m,n)] 

- ^ nr(5m|0,m)9fi [K(m,n)l 
6 

_ {3m5n) a r,, , f " ,e ,r a2 „.gm _ {SmS^} 



del [Ueirn, n)] dSi... dSm-idlWi';::^ - ^^^d^.^ [f/.(m, n)] d<5i . . . dSm-idmW^! 

J — OO J — oo 



2 



" {5l,5„)deidm [f/.(m, n)] / dSi... dSm-idmWf^. (A25) 

J — oo 

The next step is to add equations ( |A23t , l lA24b . and l |A25| l and take the continuum limit. Note that the last two lines of ( IA24t cancel with the 
last two lines of l lA25b when they are added, simplifying the end result considerably. In the continuum limit. 
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K=o(m,n) = erf - (5„0/\/2(5„ - S",,,)) 



We make use of the following properties: 



d,T^Uc=o(m,n) = —dciUc=o{m,n) = dmdciUc=o{m,n) = —djndciUc=o{m,n) 
d'^Ue=oim,n) = 9^iK=o(m,n) = -9mC'ci?7£=o(m, n) 

d,nUe=oim,n) = -9^^/6=0(^1, n). 
Rearranging the sum of equations ( IA23t , ( IA24t . and l lA25t in the continuum limit, and using the above properties yields 



dS„ P:^oiSn, S„\Sm, Sm) = -A ^'^^^^^""'"^ - B 



9™nf-(5„|0,m) 



nr(5„|o,m) 
dlnr{s^\o,m) 



nr\6^\0,m) 



dciUc=oim,n) 



(A26) 

(A27) 
(A28) 
(A29) 



(A30) 



where we have defined 



A = A{S^,S„) = {5l) - (,5^,) + 3 {5l,5„) - 3 {5^5l) 



(A31) 



B = B{S^,S„) = (S'i) + (S^Sl) - 2 {SlS„} 



(A32) 



C = C(S',„,S'„) = (5™5„) - ((5m). 



(A33) 



In orde r to evaluate equation (|A30 we need the form of the probability density nf™(5m.|0, m). To lowest order in e, this is given by equation 
(80) of iMaggiore & Riottd l l2010ij) : 



nr(5™.|0,m) 



Sm, — Sg 
TV (Sm - 5o)3/2 



exp 



(Sm — 5o)^ 



'2{Sm — So) 



Finally, substituting equation l lA34t into ( IA30b and combining with the Gaussian and Markovian term yields 



F{S„\Sm,Sm) = erk 



Sci - s„ 



+ 



\ \/'2{Sn — Sm.) 
Bi^Sm; Sn) 



cxp 



[Scl — 5mY 



1 5, 



{Scl — 5m)^ 



+ 



3S 



S71 Sm 



27r\/ Sn — Sn 



Sm 



(A34) 



(A35) 
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